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Abstract 

We consider the fermionic contributions to the free energy of noncommutative QED at finite temperature T. This 
analysis extends the main results of our previous investigation where we have considered the pure bosonic sector of the 
theory. For large values of 6T 2 (6 is the magnitude of the noncommutative parameters) the fermionic contributions 
decrease the value of the critical temperature, above which there occurs a thermodynamic instability. 
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1. Introduction 

The formulation of Quantum Electrodynamics in 
noncommutative space (NCQED) [1,2,3] opens the 
interesting possibility of new dynamics for the gauge 
fields, including self-interactions, without the intro- 
duction of additional internal degrees of freedom. 
The gauge field self-interactions of this theory con- 
tain trigonometric factors like sin^p^O^q^) where 

and q v are the momenta and O^" are the noncom- 
mutative parameters introduced via the commuta- 
tion relation of the coordinates 

[x^,x v ]*=x^*x v -x u *x fl =i6* v , (1) 

where, in general, the Gronewold-Moyal ★-product 
[1,2] between two functions f(x) and g(x) is given by 

f(x)*g(x)=f(x) expQ^>"^ g(x), (2) 

In the limit when the coordinates can be considered 
commutative, 6^ v — ► 0, all the self-interactions van- 
ish and the usual theory of non-interacting photons 
is recovered. 



When we consider a system at finite temperature 
T[4,5,6], the thermal Green functions may became 
dependent on the quantity r = 0T 2 , were 6 is the 
typical magnitude of the components 6^ v . In general 
the dependence on r may be very involved. How- 
ever, for asymptotic values of r one can have a sim- 
pler understanding of the behavior of thermal Green 
functions. It is also possible that other mass scales 
may be present, so that one could have quantities 
like 9pT , where p represents the magnitude of some 
external momenta in a given Green function. There 
is however an important quantity, namely the free- 
energy, which can only depend on r. This happens 
because the diagrams which contribute to the free- 
energy do not have external legs. 

Previous investigations on this subject [7] have 
revealed interesting properties already at the low- 
est non-trivial order (two-loop order). Corrections 
to the free energy at higher orders than two-loops, 
in thermal field theories, require for a self-consistent 
treatment to take into consideration an infinite se- 
ries of diagrams. This happens because the fields ac- 
quire, through their interactions, an effective ther- 
mal mass. This so-called plasmon effect is known 
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to occur for instance in non-abelian gauge theo- 
ries like the Yang-Mills theory at high tempera- 
ture [4,5,6]. Inasmuch NCQED is a theory of self- 
interacting gauge fields, it is natural to consider the 
possibility that similar non-perturbative effects may 
arise when we take into account all the higher or- 
der corrections to the free energy. In this case, only 
the leading hard thermal contributions are relevant. 
Consequently, other interesting phenomena such as 
the UV/IR mixing [8], which are important at T = 
0, do not arise in the present context since the tem- 
perature does provide an UV regularization of the 
thermal graphs. 

In a recent paper the free-energy of thermal NC- 
QED (TNCQED) has been analyzed by taking into 
account all the contributions which arise from the 
pure bosonic sector of TNCQED [9]. In the regime 
when r 3> 1, we have summed the ring-diagrams, 
which consist of an arbitrary number of self-energy 
insertions and yield the leading contributions. A re- 
markable property revealed by this analysis was the 
emergence of a thermodynamic instability above the 
critical temperature 



T r = 



10 « 1.229 



(3) 
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where e is the coupling constant. 

In the present paper we consider the possibility 
that the critical temperature T c may be modified by 
the inclusion of fcrmions. If we consider only Dirac 
fermions in the fundamental representation, there 
would be no modification of T c , because, as we have 
argued in [9], the contribution from fermion loops 
are the same as in the commutative theory and the 
existence of T c is a consequence of the noncommuta- 
tivc magnetic mode, which is absent in the commu- 
tative theory. However, we may additionally include 
fermions in the "adjoint" representation so that the 
fermionic action, for massless fermions, is given by 

Sfcrm = jd 4 x$*i (W - ie ltl A» * V) (4) 

where is the gauge field. The fields ip and x 
are the fundamental and "adjoint" fermions, which 
transform respectively as [3] 

V> -» e la{x) * V and x -» e ia(a ° * X * e" ia(x) , (5) 

where a(x) is the parameter of the U(l) group. In 
the limit when 9^ v — > the adjoint fermions decou- 
ple and we recover the usual electron interaction of 
commutative QED. 



In section 2 we describe the method employed 
to identify the leading contributions to the free- 
energy in terms of ring diagrams. Then, in section 
3, we present the result for the free-energy, in a form 
which allows us to derive the critical temperature, 
which turns out to be smaller than the one in the 
pure bosonic sector. The instability which occurs 
above the critical temperature may be understood, 
as discussed in section 4, in terms of the behavior of 
the noncommutative thermal magnetic mass. Many 
technical details, which are relevant in the context 
of these calculations, can be found in our previous 
paper [9]. 

2. Basic Approach 

Let us review the main steps which lead us to find 
the critical temperature T c in [9]. First we have to 
consider a non-perturbative (a very comprehensive 
analysis of the nonpertubative regime has been car- 
ried out in [10] using the lattice formulation) contri- 
bution to the free-energy which is generated by the 
ring diagrams, as follows 



&(T,0) = -- 



1 



1 ft k\ 
J VI. 



(6) 



The quantity II stands for all the contributions to 
the one-loop photon self-energy. In [9] we have shown 
that at any given order in the loop expansion, Eq. 
(6) gives the leading contribution in the limit r> 1. 
Therefore, it is implicit in (6) that we are consider- 
ing only the zero mode contribution of the photon 
propagator (denoted by wavy lines), so that k 2 = 
— \k\ 2 . The reason for this is because each photon 
propagator in (6) introduces a factor 

1 = I ( 7 ) 

k 2 (27m) 2 + |fc| 2 /V 2 ' k ' 

so that the zero mode, n — 0, yields the following 
leading contribution for large values of r 

k 2 ~ |fc| 2 { ' 

(here and in all the expressions which follows we 
have performed the rescaling (ki,k) — > T{ki,k/r) 
as well as (p4, p) — > T(p4, p) where p^ stands for the 
internal momentum of n M „) . On the other hand, pos- 
sible contributions to (6) with fermionic rings (pho- 
ton lines replaced by fermion lines) would be sub- 
leading, because the inverse of the fermionic propa- 
gator is nonzero for n — and \k\ = (the Matsub- 
ara frequencies arc all multiples of odd integers). 
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Fig. 1. One-loop diagrams which contribute to the photon 
self-energy in NCQED. 

The diagrams which contribute to IF" y are shown 
in the figure 1. The first three diagrams, involving 
only bosonic loops, have been computed in [9]. In 
the limit of large values of r, or, equivalently setting 
ki = \k\ = except inside the trigonometric factors, 
we have obtained 



2 J (2tt) 3 

OO 

^2 fWp,fc) 



where 



IV (p, k) 



p4— 2n7r 



,2 



; fcj _ kjdij, (9) 



P 



P 



(10) 



In order to compute the sum and the integral in 
(9) we have employed the following decomposition 



n 



bosc 



n boson i tt 



boson &A*fci> 



boson 



+i4 oson I v - 



k 2 ' ±A fc 2 

kuki/ k^k v 



k 2 fc 2 



(11) 



where is the heat bath four velocity and k^ = 
k a Q GVL . Substituting (9) in the left hand side of (11) 
and contracting the resulting equation with n^u" , 



k 2 J 



we can 



express the four structure functions in terms of Mat- 
subara sums and integrations over d 3 p, which can be 
done, yielding the following results for the thermal 
part 



TTboson 
ii 00 



(k) 



n 



boson 



(k) 



1 (4e) 2 T 2 

2 (2tt) 2 



1 (4e) 2 T 2 
2 



■ — coth(-7rfc) 
2k 

-coth 2 (Trk) ) , (12a) 



(27T) 2 



— coth(7rfc~ 
2k 



,i, r, 1 



__ co th>fc) + ^) (12b) 



and n^f on = n^ son = 0. 

Let us now apply the same technique to the com- 
putation of the last diagram in Fig. 1. In fact, this 
diagram represents the sum of the adjoint and the 
fundamental fcrmion loops. The contribution of the 
fundamental fermion loop is well known and gives 
the following result in the high temperature limit 

e 2 T 2 

— ir-u^Uy. (13) 



n 



fund 
pv 



As we have already anticipated this contribution of 
fermions in the fundamental representation does not 
modify the structure function II nc ; it just gives the 
usual Debye screening mass of QED [4] . The contri- 
bution of the adjoint fermion can be readily obtained 
from the Feynman rules generated by the first term 
in the action (4). Proceeding in the same way as in 
the case of the bosonic contribution, we consider the 
limit of large r, which amounts to set k± = \k\ = 
except inside the trigonometric factor. This leads to 
(reminding that the momenta have been scaled as 
P^PT) 

Ufj{k)=Ae 2 T 2 j ^1 

x E tr k^ 7 ^ (14) 




Pa = (2n+l)7r 



Using tr7 M 7 a 7^7/3 = 4 (77^^/3 + rj^-q 
we obtain 

d 3 p . 2 / k ■ p \ 



n- j (fc) = -(4e) 2 T 2 /^ )3 

OO 



sin 



p4=(2n+l)7T 



(15) 



where is given by (10). It is remarkable that 
the result for the fermionic integrand in Eq. (15) 
coincides, up to a factor, with the bosonic result in 
Eq. (10). 
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Our explicit calculation in [9] shows that there are 
two kinds of sums in terms of which the structure 
functions in (11) can be expressed, which are V?* 2 
and 1 /p 4 - The thermal part of these sums can be 
expressed in terms of the Bose-Einstein distributions 
Nb for p4 = 2mr. On the other hand, in the present 
fcrmionic case, with p4 = (2n + l)ir, the Matsubara 
sums yield a result which is simply related to the 
bosonic result in [9] by the replacement Nb — > —Np, 
where Np is the Fermi-Dirac distribution. 

Proceeding similarly to the bosonic case, we can 
now express IT^(fc) in terms of four structure func- 
tions, using the same tensor basis as in Eq. (11). 
Then, the structure functions can be explicitly com- 
puted performing the elementary angular integra- 
tions as well as the slightly more involved integra- 
tions over \p\ [11]. As a result of these straightfor- 
ward steps we obtain 



1100 _ (2tt) 2 



2 r 



TJ ad J = 



(4e) 2 T 2 
(2tt) 2 



F{k) - G{k) - 



F(k) + G(k) + 



k 2 



(16a) 



(16b) 



and = n 22 j = 0, where, for notational conve- 



nience wc have introduced 



F{ ~ k) = 7Z 

4k 




and 



G{k) = 



tanh 



irk 



- cotlr 



■nk 



(17) 



(18) 



The full result for the relevant part of the photon 
self-energy can now be obtained adding the expres- 
sions in Eqs. (16) with the corresponding ones in 
(12). Using the notation n 00 = n^ son + 0^+11^ 
and II nc = n^ son + n^J we can write the full ex- 
pression for the photon self-energy as 



n M „ = n o Uuu v + n n 



k^k v 



(19) 



3. The Free Energy 



1 T 4 



A k 



(er) 2 n o(fc) 



log 1 



Ifcl 



(er) 2 n o(fc) 



Ifcl 2 



+ 



(er) 2 n nc (fc) 



Ifcl 



(20) 



+ log 1- 



(er) 2 n nc (fc) 
Ifcl 2 



where we have introduced the quantities n o(fc) = 
n o(fc)/(eT) 2 and fi nc (fc) = n nc (fc)/(eT) 2 . This 
expression constitutes a contribution to the free- 
energy which is non-analytic in the coupling con- 
stant, since it is not proportional to a simple power 
ofe 2 . 

There are many interesting physical consequences 
which can be drawn from Eq. (20). Usually, in the 
case of commutative theories, the components of the 
static self-energy are constants independent of the 
momentum |fc|. For instance, even in the case ofther- 
mal gravity [12,13], the components of the graviton 
self-energy, in the high temperature limit, are con- 
stants (thermal masses) which can be positive or 
negative. When these constants have a positive sign, 
the free-energy picks up an imaginary part from the 
logarithm, which can be interpreted as an instability 
[14]. In the present case of the TNCQED, the self- 
energy is a function of the momentum |fc| and also 
of the direction of the vector Oi = \tijk 9jk, except 
in the case of the contribution of the fundamental 
fermion. Therefore, it is important to investigate the 
behavior of ftoo(fc) and fl nc (k) (k = |fc| sin a were a 
is the angle between fe and 6). Let us first obtain the 
asymptotic behavior for small and large k. The Tay- 
lor expansion of expressions (12) and (16) for small 
k gives 



lim n, 



no 



47T 2 - 47T 4 - 



lim n 



boson 



27T 2 r 2 87T 4 r4 
(fc) = — fc 2 - —k 4 + 



45 



945 



(21a) 



(21b) 



We have now all the ingredients to compute the 
free-energy in the limit of large r. Taking into ac- 
count the orthonormality and idempotency of the 
quantities u^u^ and and inserting Eq. (19) into 

Eq. (6) (using Y^=2 x "/ n = —x — log(l — x) in or- 
der to perform the sum of all the rings), we find 



i im n^(fc) = ^fc 2 -H^ 4 + ... . 

k^o 90 1890 



(22a) 



(22b) 
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For large k these functions behave as 



_iim fi h r n (k) = --, 

k—>oo <J 



Jim n^ son (fc) = 0, 

k^oo 



lim n$(k) 



k— >c 



2 
3' 



Jim n££(fc) = 0. 



(23a) 
(23b) 

(24a) 
(24b) 



We may consider k — ► oo as the limit when the two 
opposite charges in the dipole are very far apart 
from each other. This can be seen expressing k in 
terms of the original dimensionfull momenta so that 
k — ► |fc|r/Tsino! = 0T|fe| sin a and considering the 
limit 9T — -> oo. Both the adjoint fermion as well 
as the noncommutative photon can be pictured as 
an extended dipole of effective length L ~ 9T [15]. 
In this case, loop contributions of photons and ad- 
joint fermions reduce to the sum of two fundamen- 
tal fermion loops, corresponding to each charge in 
the dipole. Since the self-energy is proportional to 
e 2 , the two contributions add to twice the result ob- 
tained for the fundamental fermion loop. This in- 
teresting physical picture is indeed verified when we 
compare Eq. (13) with Eqs. (23a) and (24a). 

Adding all the contributions (including also the 
contribution from Eq. (13) such 
/(e 2 T 2 ) = —g), we obtain 



Tjfund 
li 00 



that ngj nd = 



lim Iloo(fc) 



Jim Iloo(fc) 



1 

3 " 
5 

~3' 



llTT 2 ,, 



47tT 4 ~ 



45 



1260 



fc 4 + 



lim n r 



90 



1890 



+ • ■ 



Jim n nc (fc) = 0. 



, (25a) 



(25b) 



(26a) 



(26b) 



We also show in the figure 2 the plots of the functions 
floo and n nc (full lines) as well as their bosonic and 
fermionic components (doted and dashed lines). 

From the behavior of n nc and n o we can now 
investigate the thermodynamic stability as follows. 





(b) 



Fig. 2. The longitudinal (a) and transverse (b) modes of the 
static photon self-energy in NCQED. 

Figure 2a shows that the contribution of the bosonic 
sector to n o does not induce an imaginary part 
to the free-energy in (20), because rioo son (doted 
line) is always negative. Similarly, the contribution 
to floo from the adjoint fermions (dashed line) is also 
negative. Hence, the longitudinal 00 mode is always 
stable. 

In the case of the transverse mode, figure 2b shows 
that the bosonic sector (doted line) is positive, yield- 
ing the instability already discussed in ref. [9]. The 
adjoint fermion contribution is also positive. There- 
fore, the sums of the fermionic and the bosonic con- 
tributions produce a larger positive result for the 
noncommutative transverse mode, which enhances 
the instability. 

Finally, we can investigate how the critical value 
of the temperature given in (3) is modified by the 
inclusion of fermions. This can be done by realizing 
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that the second logarithm in (20) has an imaginary 
part when 

( eT )2 > JL^ > JL^. (27) 

n nc (fc) " u nc (k) 

This condition can be solved analytically and leads 
to the following critical temperature 

which is associated with the noncommutative trans- 
verse mode. Comparing with (3) we conclude that 
the inclusion of the adjoint fermion has decreased 
the critical temperature. 

4. Discussion 

In conclusion, our results reveal that the possi- 
ble excitation of degrees of freedom of the adjoint 
fermion renders the system more unstable, due to 
a decrease in the value of the critical temperature, 
when compared to the one obtained in [9] . This hap- 
pens because the thermal contribution of the ad- 
joint fermion is qualitatively similar to the one of 
the (noncommutative) photon. This behavior might 
be expected since one may picture all such particles 
as having a likewise dipole structure with zero net 
electric charge [3] . Consequently, the total contribu- 
tion from the photon and the adjoint fermion to the 
transverse noncommutative magnetic mode makes 
n nc even larger, which enhances the instability of 
the system. 

One may understand this instability, if we con- 
sider the noncommutative thermal magnetic mass 
defined in a self-consistent way as 

mL = -n nc (fc = 0,fc)| fe2=mL . (29) 

Proceeding similarly as in [9] , we find that this equa- 
tion admits negative solutions for m 2 c . This occurs 
in the region er > er c , which is where the argument 
of the second logarithm in Eq. (20) becomes nega- 
tive. When er is close to its critical value, we can 
obtain an analytic solution of Eq. (29), using for n nc 
the expression given by Eq. (26a). After a simple 
calculation, we find that the noncommutative ther- 
mal magnetic mass can be written in this region in 
the form 

ml = |Ie 2 (T c 2 - T 2 ) , (30) 

where T c is given by Eq. (28). This solution shows 
explicitly that m 2 c becomes negative when T > T c . 



It is interesting to note that a negative value of the 
squared magnetic mass is reminiscent of the negative 

squared Jeans mass M 2 GT 4 which arises m 

quantum gravity [12,13]. Such a mass leads to the 
appearance of an imaginary part in the free energy, 
which indicates that the system becomes unstable 
and may undergo a phase transition [14]. 
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